One-dimensional versus two-dimensional correlation effects in the oxyhalides TiOCl 

and TiOBr 
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We have performed a comparative study of the electronic structures of the spin-Peierls systems 
TiOCl and TiOBr by means of photoemission spectroscopy and density functional calculations. 
While the overall electronic structure of these isostructural compounds is qualitatively similar, the 
bromide appears to be less one-dimensional. We present a quantitative analysis of the experimental 
dispersions in terms of exchange constant J and hopping integral t as well as a discussion of the qual- 
itative spectral features. From that we conclude that despite the one-dimensional physics triggering 
the ground state in both compounds a proper description of the normal state electronic structure has 
to take into account the anisotropic frustrated interchain interactions on the underlying triangular 
lattice. 

PACS numbers: 71.20.-b,71.27.+a,71.30.+h,79.60.-i 
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I. INTRODUCTION 

Based on magnetic measurements and LDA+U calcu- 
lations, the layered Mott insulator TiOCl was recently 
proposed to be electronically and magnetically quasi- 
one-dimensional and an unusual spin-Peierls compound. 1 
Indeed, a true spin-Peierls ground state was identified 
for TiOCl as well as for the isostructural TiOBr by 
revealing the lattice dimerization in x-ray diffraction 
measurements^^ For a conventional spin-Peierls system 
one would expect only one phase transition, however, 
TiOCl and TiOBr possess two successive phase tran- 
sitions at temperatures T C1 — 67 K, T C2 = 91 K and 
T Cl = 28 K, T C2 — 47 K, respectively. In the intermediate 
phase an incommensurate superstructure was found that 
develops below T C2 , 4,5 In a first order transition at T Cl 
this order locks in, resulting in the commensurate order 
of the spin-Peierls ground state Also in contradiction 
to a canonical spin-Peierls scenario, both compounds dis- 
play marked fluctuation effects in an extended tempera- 
ture regime above T C2 . For TiOCl, this was seen in mag- 
netic resonance? 7 ^ RamanS and infrared spectroscopy^ 
and specific heat measurements^ For TiOBr, the impor- 
tance of fluctuations was inferred from the infrared op- 
tical properties^ The origin of these fluctuations is still 
unclear but fluctuating orbital degrees of freedom can be 
ruled out due to the high energy of the lowest crystal field 
excitation in the t2 g manifold. 6-13 ' 14 From photoemis- 
sion spectroscopy on TiOCl in connection with electronic 
structure calculations we argued in a previous study that 
electronic correlations and/or spin-Peierls fluctuations 
might play an important role in this system^ The de- 
tailed behavior of the quasi-one-dimensional dispersions, 
however, is still not understood. 



In this study we compare the electronic structures of 
TiOBr and TiOCl well above the transition temperatures 
T Cl and T C2 . We present angle-integrated and angle- 
resolved photoemission (ARPES) data as well as comple- 
mentary density functional calculations. While TiOBr is 
less one-dimensional than TiOCl we conclude from a de- 
tailed analysis of the experimental dispersions as well as a 
discussion of the qualitative spectral features that proba- 
bly a Hubbard-like model including the anisotropic frus- 
trated interchain interactions on the underlying triangu- 
lar lattice is best suited to capture the relevant physics 
dominating the single-particle excitation spectra of these 
compounds. 



II. TECHNICAL DETAILS 

Single crystals of TiOCl were prepared by chemical 
vapor transport from TiCl3 and Ti02»^ TiOBr crystals 
were synthesized in a similar way from Ti02 with a 40% 
excess of Ti and TiBr 4 iii^ 

Photoemission spectroscopy (PES) was performed us- 
ing He I radiation (21.22 eV) and an Omicron EA 125 HR 
electron energy analyzer. Fresh surfaces were prepared 
by in situ crystal cleavage with adhesive tape. The re- 
sulting surfaces were clean and atomically long-range 
ordered as evidenced by x-ray induced photoemission 
(XPS) and low-energy electron diffraction (LEED), re- 
spectively. Since both TiOCl and TiOBr are insulators 
sample charging was reduced by measuring at elevated 
temperatures, i.e., at 325 K in the case of TiOCl and at 
360 K in the case of TiOBr. For TiOCl, we found from 
systematic temperature variations that the charging is 
almost negligible at and above «370K. In this situation 



2 



TiOBr 


1 1 

PES 




5 GGA+U 




A I \ 

\! | V j: 
': • Na •""*•> ; 


1 1 1 1 

TiOCI 


1 i*" 1 

A PES 




Jl/j GGA+U 






•• ... 1 .... 1 .... 1 . . 





-8 -6 -4 -2 

E-^exp (eV) 



FIG. 1: (Color online) Angle-integrated photoemission spec- 
tra and GGA+U densities of states of TiOBr and TiOCI. 

the maximum of the Ti 3d peak (see below) is located 
1.45 eV below the experimental chemical potential /i CX p, 
which corresponds to the Fermi edge position of a silver 
foil. The TiOCI spectra were aligned accordingly. The 
energy position of the Ti 3d peak in TiOBr depends on 
the temperature in a very similar way. It was aligned to 
the same value of 1.45 eV below /z eX p. The energy resolu- 
tion amounted to 60meV for the TiOCI and 110 meV for 
the TiOBr measurements. The angular acceptance was 
±1° in both cases. 

The photoemission experiments are complemented by 
calculations of the electronic structure of TiOBr and 
TiOCI within the ab initio density functional theory 
in the generalized gradient approximation with inclu- 
sion of on-site Coulomb correlations (GGA+U)i§*i£ us- 
ing the full-potential linearized augmented plane-wave 
code WIEN2k3 The GGA+U calculations, which re- 
quire a spin-polarized starting configuration, were per- 
formed by considering a ferromagnetic alignment of the 
Ti spins. For both systems a value of i?i"T max = 7 was 
used and 40k irreducible points were considered for Bril- 
louin zone integrations. The onsite Coulomb repulsion 
U and onsite exchange Jo were chosen to be 3.3 eV and 
1 eV, respectively^^ 



III. RESULTS AND DISCUSSION 

In Fig. ^ angle-integrated photoemission spectra are 
compared to the GGA+U densities of states. The calcu- 
lations yield a clear separation of the states with predom- 
inant Ti 3d character centered around — 1.5eV from the 
ones which are mainly derived from the O and Br/CI p 
levels. The latter are found between — 8.5 eV and — 3.5 eV 
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FIG. 2: (Color online) Atomically resolved GGA+U densities 
of states of (a) TiOBr, (b) TiOCI, and (c) a hypothetical 
compound with the crystal structure of TiOBr but with Br 
substituted by CI. 



for TiOBr and between -9eV and -4.5 eV for TiOCI. 
The photoemission spectra clearly show this partition as 
well. Also the significantly smaller separation of the Ti 3d 
states from the rest of the valence band which is found in 
the calculations for TiOBr is observed in photoemission. 

However, neither the exact shape nor the width of the 
Ti 3c? spectral weight are reproduced by our calculations. 
As it was pointed out in an earlier publication^ this 
is most likely due to pronounced electronic correlations 
and/or spin-Peierls fluctuations which are beyond the 
scope of the GGA+U method. A better match between 
the calculated and the experimental spectra is found with 
respect to the width of the p bands. Due to matrix el- 
ement effects a more detailed agreement of the spectra 
with the bare densities of states cannot be expected. 

In order to address the question whether the smaller 
separation between the manifold of 0/Cl(Br) and the Ti 
3d states in the Br compared to the CI system is caused 
by the change in chemistry or rather related to the struc- 
tural expansion we performed a third calculation. There 
we assumed a hypothetical compound with the crystal 
structure of TiOBr but with Br substituted by CI. In 
Fig. |21 the atomically resolved densities of states for all 
three systems is presented. The decomposition into Ti 
states near the chemical potential on the one hand and 
O and Cl(Br) states at higher binding energies on the 
other becomes very clear in this plot. By comparison 
it is easily seen that the separation of the Ti states is 



FIG. 3: (Color online) EDCs of TiOBr along the crystal- 
lographic 6 axis, corresponding to the YrY line in the or- 
thorhombic Brillouin zone. 



smaller in the Br system mainly due to the expanded 
unit cell of TiOBr. 

Turning to our angle-resolved photoemission data, we 
display energy distribution curves (EDCs) of TiOBr in a 
broad energy range along the crystallographic b direction 
in Fig. 01 What immediately stands out are the well- 
pronounced dispersions, especially in the Br Ap/O 2p part 
of the spectra. The dispersions are clearly symmetric 
with respect to the T point reflecting atomic long range 
order and thus can be taken as a further proof of the 
good surface quality. 

In Fig. H the Ti 3d spectral weight distributions of 
TiOBr and TiOCl are compared. The four panels show 
ARPES intensity plots I(k, E) along the crystallographic 
axes a and b. As seen, small but clear dispersions are 
discernible in all four panels. These are quantitatively 
identified by lines and markers. The lines indicate EDC 
peak maxima obtained by a fitting procedure. For both 
compounds, this works well along the a direction and in 
the first half of the path from T to Y, i. e. along b in 
direct space. In the second half, however, the EDCs lose 
their single peak-like shape as an additional feature at 
higher binding energies appears and thereby the spectra 
gradually change their form to end up in a broad hump 
at the Y point. 13 Thus a fit with a single peak is not 
adequate here. Instead, we take the first order moment 
of the EDCs at the Y point, indicated by crosses in Fig. 01 

Following first the dispersions along the a axes, i. e. 
from r to X in reciprocal space, both compounds dis- 
play a qualitatively similar behavior: Initially the en- 
ergy changes only slightly before it reaches a maximum, 
and then shifts downwards until the Brillouin zone edge. 
Also along the crystallographic b direction, the disper- 
sions show a qualitative resemblance in the two com- 
pounds. Starting from T, a single peak disperses up- 
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FIG. 4: (Color online) ARPES intensity plots I(k,E) of 
TiOBr and TiOCl along the crystallographic axes a and b, 
corresponding to the XFX and YFY lines in the orthorhom- 
bic Brillouin zone. Blue lines indicate the peak positions of 
the EDCs. The blue crosses mark the first moment of the 
EDCs at the Y points. 



TABLE I: Dispersions w a , Wb, and wy of the Ti 3d weight 
in TiOBr and TiOCl, measured from the maximum to the 
minimum energy of the peak for the direction a, and from the 
maximum energy of the peak to the first moment of the Y 
point EDC for the direction b, respectively. 



width 


TiOBr 


TiOCl 


W a 


0.27(3) eV 


0.12(3) eV 


Wb 


0.26(5) eV 


0.47(5) eV 


w'b 


0.13(1) eV 


0.17(1) eV 



wards, with the shape of the EDCs remaining essen- 
tially unaltered until ^rY. As already pointed out, sub- 
sequently the spectral changes get more complex. In any 
case the overall dispersion bends downwards reaching a 
minimum at the Y point which lies further below the one 
at r. We emphasize once more that while at first glance 
the dispersions from the intensity plots might appear to 
be similar along the a and b directions they are qualita- 
tively different if one looks at the spectral shapes of the 
corresponding EDCs. 

For a quantitative comparison of the electronic disper- 
sions between the two compounds we list the parameters 
W a , Wb, and w' b in TableQ] The unprimed quantities refer 
to the overall dispersion widths while w' b measures the 
width of only the inner part of the dispersion in the re- 
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gion from T to about ^TY. These widths are determined 
either from the difference of the maximal and minimal 
peak energies along the corresponding paths in k space 
as obtained by the fitting procedure (see lines in Fig.0J, 
or the difference between the maximum peak energy and 
the first order moment at the Y point in case of Wb (see 
markers in Fig.^J. The errors indicated reflect the scat- 
ter from several samples and measurements. 

As is immediately read off from Table [I] the overall dis- 
persion width along b is significantly smaller in TiOBr 
with respect to TiOCl while it is the other way round 
regarding the a axis. Given that both systems and in 
particular their electronic structures are governed by the 
same physics we hence conclude that the anisotropy of 
TiOBr is less pronounced than in TiOCl. This is in line 
with the trend in the relevant hopping integrals as de- 
rived from downfolding LDA+U results as well as with 
the fact that the Bonner-Fisher curve indicative for ID 
Heisenberg chains does not provide a good fit for the 
high-temperature magnetic susceptibility of TiOBr in 
contrast to TiOCl32i 

At this stage we resort to the central result of our previ- 
ous study on TiOCli^ where among various approaches to 
the k resolved electronic structure we identified the spec- 
tral function of the single-band ID Hubbard model as 
the most promising starting point in order to get further 
insight into the microscopic physics behind our experi- 
mental data. We recall that the most obvious shortcom- 
ings of the model calculations were the lack of any evi- 
dence for spin-charge separation and the absence of the 
so-called shadow band in the ARPES spectra although 
for the relevant parameter regime it should be clearly 
observable^ As possible explanations for the discrepan- 
cies between theory and experiment wc invoked multi- 
orbital effects and/or spin-Peierls fluctuations. A com- 
parison between the isostructural compounds TiOBr and 
TiOCl now opens the possibility to further explore the 
reasons for this discrepancy between ID Hubbard model 
prediction and experiment. For this purpose we focus on 
the central part of the dispersion along the b axis marked 
by lines in Fig. 01 Within the single-band ID Hubbard 
model this part corresponds to the (u>, k) region of the 
spinon and holon branches (see Fig. 7 in Ref . fl3f) . Hence, 
depending on the dominant character of the experimen- 
tal dispersion its width w' b should either scale with the 
exchange constant J whose value can be extracted from 
magnetic susceptibility measurements^ or the hopping 
integral t. The latter either can be inferred from the 
Hubbard model perturbation expression for the exchange 
constant, J = At 2 /U , and thus should scale as t oc \fj 
or it can be deduced from an appropriate downfolding 
procedure of LDA+U band calculations. 20 An account 
of these quantities is given in Table |Hj The values of 
J for TiOBr and TiOCl differ by -45% (with respect to 
the value in TiOCl), whereas the experimental dispersion 
width is smaller in TiOBr by only 23%. 

On the contrary, the experimental width w' h nicely 
matches the transfer integral as obtained from t = 



TABLE II: Comparison of TiOBr and TiOCl with respect to 
various parameters possibly relevant for the observed disper- 
sions: Width of the central part of the Ti 3d spectral weight 
dispersion along b (w' b ), exchange constant J obtained from 
the magnetic susceptibility,- hopping probability t derived 
from J, and the hopping integral t of downfolding LDA+U 
studies p— 





TiOBr 


TiOCl 


Ax/x 


ID dispersion w' b 


0.13 eV 


0.17 eV 


-23% 


J from magnetic susceptibility 


32meV 


58 meV 


-45% 


t = V7U/2 


0.16 eV 


0.22 eV 


-26% 


t from LDA+U 


0.17 eV 


0.21 eV 


-19% 



V JU /2 which is smaller in the bromide by 26% compared 
to the chloride. A similarly fair agreement is achieved 
with t from the LDA+U downfolding studies 2 ^ where the 
effective d xy -d xy hopping parameter along the b axis is by 
19% smaller in TiOBr with respect to TiOCl. From this 
analysis it follows that the experimentally observed dis- 
persions scale with t — not with J — and thus can clearly 
be identified as charge excitations. Moreover, there is no 
indication in the spectra of Fig. [3| for any asymmetry to- 
wards lower binding energies which could be interpreted 
as a remnant of the spinon branch. Thus, we are led to 
rule out that the experimentally observed spectral dis- 
persion is a superposition of both holon and spinon exci- 
tations. 

Turning back to our starting point, i.e., the single- 
band ID Hubbard model, a simple explanation for the 
lack of a spinon branch and shadow band would be that 
the two-dimensional coupling in these compounds under 
the particularities of a triangular lattice is already large 
enough so that generic ID features of the spectral weight 
distribution can not persist. One would then be left with 
a lower Hubbard band which though incoherent in na- 
ture could still display sizable dispersion alike to what 
we observe experimentally. Clearly, this issue demands 
for more detailed theoretical investigation. Alternatively, 
one could stick to the ID Hubbard model and cite the 
coupling to phonons, multiorbital or spin-Peierls fluctu- 
ation effects as possible causes for the complete suppres- 
sion of the generic ID phenomenology. However, from the 
organic quasi-one-dimcnsional conductor TTF-TCNQ we 
know that the phonons which should couple equally 
strong to the electrons in this charge density wave system 
and are of even higher energy, do not suppress or smear 
out the spectral features of spin-charge separation and 
shadow band. Moreover, since as already stated the or- 
bital degrees of freedom are quenched in the oxyhalides 
and spin-Peierls fluctuations should not be too impor- 
tant so far above the transition temperatures, these two 
effects are probably not very effective. 

In view of the above results we are thus left in a sit- 
uation where on the one hand the spin-Peierls ground 
state of both oxyhalides studied here is dominated by 
ID interactions while the analysis and discussion of the 



5 



electronic dispersions at room temperature as well as the 
incommensurate order in the intermediate phase point 
to the importance of 2D (frustrated) interchain interac- 
tions. It remains interesting to see whether or not an 
anisotropic Hubbard model type description taking ac- 
count of the magnetic interchain frustrations on the un- 
derlying triangular lattice is capable of better describing 
the electronic properties of this class of materials posi- 
tioned in the regime between ID and 2D correlations. 

IV. CONCLUSIONS 

In summary, we have investigated the electronic struc- 
ture of TiOBr and TiOCl by means of photoelectron 
spectroscopy and GGA+U calculations. In these com- 
pounds we find a very similar density of states, with the 
exception of an enlarged separation between the manifold 
of O and Cl(Br) levels and the Ti 3d states. This differ- 
ence is due to the enlarged unit cell in TiOBr. ARPES 
measurements reveal that the Ti 3d spectral weight dis- 
persions of the two compounds show a great qualitative 



resemblance, with the main difference that in TiOBr the 
overall dispersion width is larger along the a axis and 
smaller along the b axis, compared to the chloride. Thus 
the quasi-one-dimensional character of TiOBr is clearly 
less pronounced. From the comparison of TiOBr and 
TiOCl, we are able to show that the dispersion of the 
d band along the b axis scales with the square root of 
the exchange constant J, i.e. with the hopping integral 
t. Based on this analysis and the lack of any sign for 
spin-charge separation or shadow band and against the 
background of other experimental findings we argue that 
the physics of both systems is governed by the interplay of 
both one-dimensional and two-dimensional correlations. 
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